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The doubling of the superorbital period of Cyg X-l 
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ABSTRACT 

We study properties of the superorbital modulation of the X-ray emission of Cyg X-l. We find 
that it has had a stable period of ~300 d in soft and hard X-rays and in radio since 2005 until at 
least 2010, which is about double the previously seen period. This new period, seen in the hard 
spectral state only, is detected not only in the light curves but also in soft X-ray hardness ratios 
and in the amplitude of the orbital modulation. On the other hand, the spectral slope in hard 
X-rays, > 20 keV, averaged over superorbital bins is constant, and the soft and hard X-rays 
and the radio emission change in phase. This shows that the superorbital variability consists 
of changing the normalization of an intrinsic spectrum of a constant shape and of changes of 
the absorbing column density with the phase. The maximum column density is achieved at the 
superorbital minimum. The amplitude changes are likely to be caused by a changing viewing 
angle of an anisotropic emitter, most likely a precessing accretion disc. The constant shape 
of the intrinsic spectrum shows that this modulation is not caused by a changing accretion 
rate. The modulated absorbing column density shows the presence of a bulge around the disc 
centre, as proposed previously. We also find the change of the superorbital period from ~150 
d to ~300 d to be associated with almost unchanged average X-ray fluxes, making the period 
change difficult to explain in the framework of disc-irradiation models. Finally, we find no 
correlation of the X-ray and radio properties with the reported detections in the GeV and TeV 
y-ray range. 

Key words: accretion, accretion discs - radio continuum: stars - stars: individual: Cyg X-l 
- stars: individual: HDE 226868 - X-rays: binaries - X-rays: stars. 



1 INTRODUCTION 

A major characteristic time scale in X-ray binaries is their or- 
bital period, which, in many cases, is seen in their l ight curves 
as periodic modulations, see, e.g., IWen et alj (2006) for a list 
of orbital periodicities from the All- Sky Monitor (ASM) aboard 
Rossi X-ray Timing Explo rer (RXTE; Bradt, Rothschild & Swank! 
1 19931 : lLevineetalj|l996l) 



pear 



due to a number 



The observed flux modulations ap- 
of differen t physical effects, see 
IPoutanen. Zdziarski & Ibragimovl d2008h (hereafter PZI08) for a 
list of the possibilities. 

In addition, a number of X-ray binaries show modulation 
at quasi-periods much longer than their orbital per i ods. Those 
long-period modulation s, listed, e.g., in IWen et all j2006h and 
Ogilvi e & Dubusl fcOOlh . hereafter OD01, are called superorbital 
(hereafter SO). In most of the known cases, the SO periodicity 
(or quasi-periodicity) appears to be caused by precession of an ac- 
cretion disc and/or jet, which either results in variable obscuration 
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of emitted X-rays as in Her X-l jKatzlll973l : iRobertsll 1974) . or in 
changes the viewing an gle of the presume d anisotropic emitter, as 
in SS 433 <Katdll98Ch or Cyg X-l (e.g.. iLachowkz et alj|200d 
hereafter L06. Ilbragimov. Zdziarski & Poutanenl 120071 hereafter 
IZP07), or both. In some relatively rare cases, the SO periodicity 
is caused by modulation of th e accretion rate, e.g., in the low-ma ss 
X-ray binaries 4U 18 20-303 Izdziarski Wen & Gierlinskill2007 
and 4U 1636-536 dShih et alj 120051 : iFarrell, Barret & Skinne 
l2009h , 

In particular, Cyg X-l, a well-known high-mass X-ray bi- 
nary containing a pers istent accreting black-hole source (see, 
e.g., [ziolkowski] 120051 for the properties of the binary), has 
shown such (quasi) periodicity in its X-ray, optical and ra- 
dio emission with an average SO period o f P sm - 150 
d, found in a large number of investiga t ions dBrocksopp et alj 
1 999 a:; IPoolev. Fender & Brocksopd 1 1 999 ; Kitamoto et al. 2000; 



Karitskava et alj|200ll ; l6zdemir & Demircarj|200ll ; Benlloch et alj 
200 ll l2004l :~L06; IZP07). The m odulation is detectable in the 



hard X-ray spectral state (see, e.g., [Zdziarski & Gierlinskill2004|) . 
in which Cyg X-l spends most of the time. No SO modulation has 
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Figure 1. (a) The light curves starting from MJD 53500 showing the daily averages from ASM (black, 1.5-12 keV), BAT (red, 15-50 keV), and the Ryle/AMI 
radio telescope (green, 15 GHz), normalized to their respective average values within the hard state of MJD 53880-55375 (delineated by the vertical dashed 
lines). For clarity, the uncertainty of the radio points is not shown. The vertical dotted lines delineate the previous hard state, MJD 53540-53800. The arrows 
show the times of the reported y-ray detections by MAGIC (magenta) and AGILE (black), (b) The corresponding 3-12 keV photon spectral indices, determined 
using the ASM data. The data at MJD > 55375 show the 2010 soft state of Cyg X-l. 



been detected in the soft states, but this seems to be due to their 
typical short duration and strong overall variability. The underly- 
ing physical process appears to still operate, as indicated by the 
reappearance of the SO modulation at a phase expected from past 
modulation at a next hard state (L06; IZP07). 

The reality of the ~150-d period was then questioned bv lRicol 

(2008), hereafter R08. Although he did find the strongest, and 
highly significant, power-spectrum peak at 148 ± 2 d using data 
similar to those analysed by IZP07 (who found P sup ^152 d), he 



claimed it to be "an artefact of Fourier-power-based analysis". The 
basis for this claim was his finding of a variable P sllp in three sep- 
arate sub-samples each covering ^300 d. We disagree with this re- 
jection of the reality of the 150-d peak. First, determination of the 
modulation period based on a single interval covering = 2P sup is 
obviously questionable. Although R08 claimed those individual pe- 
riods were highly significant statistically, he determined their prob- 
ability under the assumption of the presence of white noise only, 
which is well known to strongly overestimate the actual signifi- 
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cance in astrophysical sources (e.g jlsrael &StellJ 1996b . Second, 
unlike the periodic orbital modulation, the SO modulation repre- 
sents a quasi-periodic modulation, in which the underlying clock 
is not perfect, and where the period changes from one cycle to an- 
other. Furthermore, the presence of an aperiodic long-term variabil- 
ity, related to the variable stellar wind from the compan ion (which 
also causes spectral state changes in Cyg X-l, e.g., iGies et all 
l2003h . and possible accretion disc instabilities, obviously distorts 
the SO period (presumably related to precession) if measured in a 
short interval. 

Thus, the 150-d period does represent the average P sup as mea- 
sured from ~ 1980 to 2004 (L06; IZP07). The physical reality of 
this average period was then strongly confirmed by PZI08, who 
found that both the X-ray spectral hardness and the depth of the 
X-ray orbital modulation are modulated with this P sup . 

On the other hand, R08 found a striking change of the SO 
period to P sup >300 d occurring in 2005. The data at his dis- 
posal, until 2008, covered only ~ 3-P sup , during which the pe- 
riod was = 320-330 d. Notably, the new period is similar to that, 
294 ± 4, seen in the optical and X - ray ranges from 19 69 to ~1980 
jPriedhorskv, Terrell & Hoitll 1983k Kemp et "allll983h , and appar- 
ently in the optical data till 1985 fempet all 19871) . 

In the present work, we find the ~300-d period t o persist un- 
til the recent switch to th e soft state in 2010 June l lNegoro et al.l 
l20ld : lRushton et al.ll20ldh . We find the new P sup ~300-d modula- 
tion to have the same physical properties as those found by PZI08. 
Namely, it appears not only in the soft and hard X-ray and radio 
light curves, but also in the phase dependence of the X-ray hardness 
and of the depth of the X-ray orbital modulation. Furthermore, we 
find very strong evidence that the SO modulation is due to changes 
of the source geometry but not due a variable accretion rate. 



2 THE DATA 

We study monitoring data from the RXTE/ 'ASM, which contain en- 
ergy channels at 1.5-3 keV, 3-5 keV and 5-12 keV. We apply a 
small correction for the empty-field ASM count rates of 1 mCrab, 
appropriate for Cyg X-l (R. Remillard, personal communication). 
We then use the mon i toring by the Burst Al ert Telescope (BAT, 
iBarthelmv et "all l2005l ; iMarkwardt et alj|2005h on board the Swift 
satellite. We use both the 15-50 keV data from the BAT web pagqj 
and 14-195 keV 8-channel light curves created for this work. The 
channels are between energies of 14, 20, 24, 35, 50, 75, 100, 150 
and 195 keV. 

The ASM c ount rates are convert ed into energy fluxes us- 
ing the matrix of Zdziarski et al.l d2002l) . hereafter Z02. We have 
checked that the fluxes obtained are within 10 per cent of those ob- 
tained by scaling the count rates to those from the Crab and using 
the best fit to the 0. 5-180 keV Crab energy spectrum to data from 
several detectors of iRao & Vadawalel ( |200lt) . We convert the BAT 
8-channel data into energy fluxes by scaling to the Crab spectrum. 
We use the results of lRao & Vada wale ( 2001 ]) up to 100 keV , and at 
higher energies we use those of Ijourdain & Rogu es (2009), based 
on the INTEGRAL/SPl data. This gives us, 



F(E) = A, exp(-/V H cr)(£/l keV) 



-r,+i 



(1) 




where the H column density is /V H = 3.23 ± 0.02 cm 2 , <r is the 



http://swift.gsfc.nasa.gov/docs/swift/results/transients 



Figure 2. The Lomb-Scargle periodograms for the ASM 5-12 keV dwell 
data (red curve) obtained using the hard-state intervals of MJD 53540- 
53800 and 53880-55375. The horizontal dashed line shows the level of the 
90 per cent significance for . We see many apparently highly significant 
peaks. The black curve shows the ASM periodogram divi ded by the one 
descri bing the red noise (blue curve) using the method of Israel & Stellal 
1 19961) . see Sectionfjl] Now only the -300 d and the (orbital) 5.6-d peaks 
are significant. 



bound-free cross section per H atom averaged over the cosmic com- 
position, A, = 10 ± 1 cm" 2 s~\ Ti = 2.10 ± 0.01 for E < E h = 100 
keV, Ti = 2.24±0.02 at E > E\y, and Ao follows from the continuity 
at£ b . 

In addition, we use 15-GHz data from the Ryle Telescope and 
the AMI Large Array. The AMI Large Array is the re-built and 
reconfigured Ryle Telescope. The monitoring by Ryle was carried 
out until 2006 June 16, and the AMI monitoring of Cyg X-l started 
2008 April 17. The new correlator has a bandwidth of about 4 GHz 
(compared with 350 MHz); but the effective centre frequency is 
similar, and in any case the radio spectru m of Cyg X-l, at le ast in 
the hard state, is known to be very flat dFender et al.ll200d) . The 
data are subject of variations in the flux calibration of about 10 per 
cent from one day to another. 

We take into account the barycentric correction for all the data. 
We first investigate data from the time interval considered by R08 
but extended until 2010. Also, we study the uninterrupted long hard 
state of 2006-2010, which lasted -1500 d. 



3 RESULTS 

3.1 The light curves and power spectra 

Fig.QI a ) presents the daily-average light curves from ASM (1.5- 
12 keV), the BAT transient monitoring (15-50 keV), and the 
Ryle/AMI 15-GHz monitoring. They have been renormalized to 
their respective average values during the long 2006-2010 hard 
state, MJD 53880-55375. The (un-weighted) averages are 20.6 s _1 , 
0.173 s" 1 , and 12.7 mjy, respectively. The corresponding rms, i.e., 
the fractional intrinsic standard deviation (after subtraction of the 
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Figure 3. (a) The MJD 53880-55375 light curves from ASM at 1.5-3 keV (black Mangles), 5-12 keV (blue circles) and BAT at 15-50 keV (red crosses) 
folded over the SO period for the ephemeris of equation (3)- The dependencies are normalized to their respective average values, (b) The corresponding photon 
spectral indices in the energy ranges of 1.5-5 keV (black triangles), 20-35 keV (red crosses), 35-75 keV (magenta squares) and 75-150 keV (blue circles), (c, 
d) The light curves and spectral indices as in (a, b), respectively, but folded over the orbital period. 



contribution from the measurement errors), is 0.27, 0.19 and 0.35, 
respectively. 

We see that the fluxes from the two X-ray instruments closely 
follow each other during the 2006-2010 hard state as well as the 
previous one, MJD 53540-53800. The radio flux, when available, 
strongly correlates with the ASM flux in the hard st ate. (We discuss 
the ra dio/X-ray correlation in Cyg X-l in detail in lZdziarski et al.l 
1201 ll .) We also calculate an equivalent spectral index for each pair 
of adjacent energy channels using the method of Z02, see Appendix 
|A]for details. This spectral index corresponds to a power law spec- 
trum that would have the same energy flux in each of the channels 
as that of the data. During both hard states, the 3-12 keV photon 
spectral index was T ~ 1.7, and it was always (except for three 
isolated occurrences) < 2.1, as shown in Fig.[]2b). We see a signifi- 
cant increase of the ASM flux around MJD ~55000, during which, 
however, the source remained in the hard state, as shown by T < 2 
and the correlation of the ASM and the radio fluxes. 

Fig. [T] also shows (arrows) the times when the Major At- 
mospheric G a mma Imaging Cherenk ov Telescope ( MAGI C) 
dAlbert et alj 120071) and AGILE dSabatini et al] l2010j|bt 
Bulgarelli et al. 2010) reported detections of Cyg X-l in TeV 
and GeV respectively. We see that these events do not seem to 
correlate with any source properties in X-rays, except for the last 
detection by AGILE, which took place around the transition to the 
soft state in 2010 June. All of the AGILE detections occurred in 
rather low radio states. 

In Fig.[T] we also clearly see the ~300-d modulation (exclud- 
ing the soft states), covering now 6 periods. The significa nce of 
this pe riodi city can be qua ntified using the periodogram of iLombl 
i 19761) and lScargld dl982h . The periodogram obtained using the 
ASM 5-12 keV (which channel is least affected by the orbital mod- 
ulation) dwell data is shown by the red curve in Fig. [2] The peri- 
odogram for the BAT 15-50 keV orbital data is similar, and, for 
brevity, it is not shown. Under the assumption of the presence of 



periodicities and white noise only, the probability of finding a peri- 
odicity by chance in one of N trials is, 

p N = l-(l- Pl ) N , (2) 

where p\ = exp(-PLs) is the probability for a single trial, Pls is the 
Lomb-Scargle power and the number of independent frequencies, 
N, approximately equals 1/2 of the number of the independent data 
points. Using this definition, we would have obtained many peri- 
odicities with extreme significance. This was, in fact, done by R08, 
who claimed, e.g., the chance probability of 10~ 96 for the ~300-d 
peak in his ASM periodogram (obtained using daily averages, and 
until MJD 54592). 

However, it is well-known that power spectra of X-ray bina- 
ries, and in particular, Cyg X-l, show quite strong components at 
low freq uencies, and are th us not at all compatible with white noise, 
see, e.g. jReig et al. I J2002h . This effect, often called red noise, can 
be taken into account, e.g., by averaging the periodogram over fre- 
quencies surrounding a given peak, and th en dividing the origi - 
nal Lomb-Scargle power by the average one (Israel & Stella 1996). 
We use here this method as implemented in 1 Wen et al.l ( 120061) , see 
that paper for details. The corrected probabilities are now given by 
the same formula as above, but using the rescaled P L $, which is 
shown by the black curve in Fig. [2] For the ASM and BAT data, 
the rescaled peak power has the values of .P LS =i 20, 16.5, and is 
at ^sup - 319 ± 5 d, 312 + 5 d, resp ectively, where the uncertainty 
has been calculated using eq. (3) of IWen et al.l d2006h . For the re- 
spective N = 15174 and 10462, we find p N = 3 x 10~ 5 , 7 x 10~ 4 . 
We note that these are conservative estimates since the number of 
independent data points is significantly lower (and the correspond- 
ing probability lower) because the aperiodic components of light 
curves clearly exhibit a degree of autocorrelation, and thus adja- 
cent measurements within the autocorrelation length are not inde- 
pendent, see, e.g., L06. However, since we have already established 
the significance of the periodicity at a very high level, we ignore 
this complication. Note that we do not find any peak at P sup = 1000 
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d, which was found by R08 and speculated to be real, but which 
appears to simply reflect the total duration of his data set. 

Given the 2006-2008 gap in the radio coverage, the power 
spectrum of the radio data has a rather broad SO peak, and we do 
not study it in detail. Still, the maximum LS power corresponds to 
Psup = 316, in agreement with those of the ASM and BAT data. 

The full data set considered above suffers from being inter- 
rupted by the occurrence of the soft state during MJD 53800- 
53880, which appears to have introduce some phase shift into the 
SO modulation resulting in an apparent lengthening of the first cy- 
cle, see Fig.[TJ On the other hand, the data sets for MJD 53880- 
55375 offer us an opportunity to study a very long, uninterrupted, 
hard state covered by both the ASM and BAT as well as partly 
by the radio monitoring. The periodograms for MJD 53880-55375 
look similar to those shown in Fig.[2] and thus we do not show them 
here. We find P 6up = 315 ± 6 d, 306 ± 6 d for the ASM 5-12 keV 
and BAT 15-50 keV data, respectively. 
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Figure 4. The distribution of the fractional frequency of occurrences of X- 
ray dips (see Section l3~2l over the orbital phase during MJD 53880-55375 
(black squares) and the entire currently available ASM data, MJD 50087- 
55406 (blue triangles). The uncertainties shown assume Poisson's statistics. 



3.2 Folded light curves 

For the MJD 53880-55375 data set, we now search for the de- 
pendence of the X-ray hardness and the depth of the orbital X-ray 
modulation on the SO phase, which effects were earlier found by 
PZI08 for the ~150-d SO modulation. We average now the loga- 
rithms of the count rates (without weights), which is motivated by 
their normal distributions, whereas the count rates themselves have 
the distributions far from normal (PZI08). The uncertainties on the 
averages are calculated from the dispersion of the individual mea- 
surements (without taking into account their errors). As above, we 
use the ASM dwell and BAT orbital data. The zero SO phase, O = 
is taken to correspond to the minimum ASM count rate. To facil- 
itate the study of the dependence of the orbital modulation on the 
SO phase, we take a value of the SO period within the uncertainties 
of the periods found in the data but equal to a multiple of the orbital 
period, P orb = 5.6 d. This results in the ephemeris, 

min[MJD] = 53275.8 + 313.6£, (3) 

where E is an integer. The resulting folded curves of the count 
rates are shown in Fig. [3f a )- We see that the amplitude of the SO 
modulation decreases with energy. The amplitude at 1.5-3 keV 
is 0.6 < FJ(F) < 1.3 whereas it is only about ±0.1 in the 15-50 
keV range. This is compatible with the modulation being enhanced 
by bound-free absorption around the flux minima (as found by 
PZI08 for the ~150-d modulation), as well as with the degree of 
anisotropy of Compton-scattered spectra decreasing with photon 
energy, see IZP07. Fig.[3jb) shows the corresponding spectral hard- 
ness as given by the photon spectral index calculated using the 
method of Z02. We see that indeed the 1.5-5 keV spectra become 
significantly harder at the minimum of the SO modulation. On the 
other hand, we see no spectral modulation above 20 keV, as shown 
by three spectral indices in the 20-150 keV range. The SO modula- 
tion thus consists of a modulation of the total flux and of the column 
density only, but not of the slope of the intrinsic spectrum. This is 
fully compatible with the modulation being caused by changing the 
viewing angle of the precessing accretion disc, but not by changing 
the accretion rate. 

Compared to the previous hard states covered by the ASM, 
from MJD 50350 to 53 174 (IZP07; PZI08, during which period was 
^sup - 150 d), the present SO modulation is significantly stronger, 
especially the minima caused by absorption are deeper. The SO 
modulation of T( 1.5-5) keV is also much stronger now, as com- 



pared to its variability within T ^ (1.3-1.5) during the previous 
hard states. 

On the other hand, we find the ASM count rates averaged over 
the hard states characterized by P sup a 150 d and P sup ^ 300 d to 
be to be almost identical. The former in the 1.5-3 keV, 3-5 keV 
and 5-12 keV bands are higher than the latter by only 7, 2 and 3 
per cent, respectively. The largest difference being in the 1.5-3 keV 
rate is compatible with the found increased absorption during the 
P sup - 300-d epoch, whereas the rates at to higher energies (corre- 
sponding to the intrinsic emission and not significantly affected by 
absorption) are practically identical. Furthermore, the bolometric 
fluxes calculated using the ASM and Compton Gamma Ray Obser- 
vatory BATSE data, and the ASM and BAT data, which correspond 
to these two respective e pochs, are also almost identical, as found 
bv lZdziarski et alj J20TTI) . 

Figs.[3lc,d) show the same light curves and indices as in (a,b) 
b ut folded over the 5 .6 d orbital period, fo r which the ephemeris 
of lLaSala et all Jl998h and lBrocksopp et al.1 dl999bl) has been used 
(see eq. [1] in IZP07). Overall, the orbital modulations are quite 
similar to the SO ones. The main difference is the relative weakness 
of the 15-50 keV modulation (though still significant), now about 
±2 per cent. This modulation is also found in the periodogram, 
which has the rescaled Pls - 6.5 at P or b = 5.6 d, for which the 
single-trial probability (appropriate for modulation at a known fre- 
quency, see, e.g., L06) is p\ =; 0.0015. These results are similar 
(though more significant) to those found by L06 for the BATSE 
20-100 keV data. The modulation is explained (L06) as being due 
to Compton scattering of hard X-rays away from the line of sight 
by the same gas (of the stellar wind of the donor star) as that caus- 
ing th e orbital m odulation of soft X-rays due to bound-free absorp- 
tion dWen et alj fl999). Again, we see no modulation of the spec- 
tral indices at E > 20 keV, which is explained by the weakness 
of the bound-free absorption and constancy of the Compton cross 
section at those energies. Compared to the previous hard states 
(MJD 50350-53 1 74), the present orbital modulation is significantly 
stronger, e.g., the minimum of the 1.5-3 keV relative flux was be- 
fore about 0.85, and r(1.5-5 keV) was within 1.3-1.6 only. 

Another measure of the orbital X-ray modulation is the fre- 
qu ency of X-ray dips, defined by i ncreases of the X-ray hardness 
by iBalucinska-Church et a They are presumed to origi- 

nate dur ing passage of clouds of the clum py wind through the line 
of sight. Baluciriska-C hurch et alj J2000h defined them as events in 
the ASM dwell light curve with the ratio of counts of either (3- 
5 keV)/(1.5-3 keV) > 2 or (5-12 keV)/(3-5 keV) > 2.5. Fig. H 
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Figure 5. (a) The MJD 53880-55375 light curves from the ASM 1.5-3 
keV folded over the orbital period and for a chosen interval of the SO phase 
of equation (3): <D = -0.05 to +0.05 (black triangles), and 0.55-0.65 (red 
squares). The dependencies are normalized to their respective average val- 
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Figure 6. (a) The MJD 53540-53800 and 53880-55375 light curves from 
the Ryle/AMI radio telescope folded over the SO period of P sup = 313.6 d. 



shows the fractional occurrence of the dips as a function of the 
orbital phase. We see the dips are concentrated around th e orbital 
phase 0, as found before dBalucinska-Church et al ] |200d : PZI08), 
with some shift of the peak towards > 0. We also found an SO 
dip modulation, with similar shape as in PZI08. Compared to the 
previous hard states, within MJD 50087-53789, studied by PZI08, 
we find the dips are now much more frequent. We find 2251 of 
them among 25147 of the data points, compared to 1151/31211 
in the study of PZI08. In Fig. [4] we also show the dip distribu- 
tion for the entire currently available ASM data, with the frequency 
of 3929/96018 (similar to that of PZI08). The present increase of 
the dip frequency is compatible with the orbital modulation being 
stronger during the presently considered interval than at previous 
occurrences of the hard state (see above). In either case, we see 
no evidence of the secondary pe ak at </> - 0.6 claimed (based o n 
data with rather low statistics) by Baluciriska-Church et a[| |2000), 
in agreement with the finding of PZI08. 



PZI08 found that the orbital modulation is significantly 
stronger near the minimum of the SO cycle than near its maxi- 
mum, and fitted its dependence using a model of a bulge around 
the centre of the precessing disc. Here, we confirm it for the present 
data (with the higher P sllp ). Fig.[5ja) shows the profile of the orbital 
modulation for two ranges of the SO phase; around the SO mini- 
mum, O = (-0.05 to +0.05), and near the maximum (see Fig.(3jt), 
O = (0.55-0.65). We see that the latter profile is much shallower, 
0.85 < FJ{F) < 1.1, than the former, +0.4. This is also very clearly 
seen in the modulation of the 1.5-5 keV spectral index, see Fig. 
Ob), again much stronger at the former range of O than the latter. 

Fig. |6] shows the SO modulation of the radio flux. Given the 
2006-2008 gap in the radio coverage, we use now both of the occur- 
rences of the hard state. We see a very significant modulation, again 
stronger than during the previous hard states (which relative ampli- 
tude was about 0.7-1.3). As discussed in IZP07, the most likely 
cause of the modulation is the precession of the jet tied to that of 
the disc. Then, the jet emits anisotropically due to the relativistic 
beaming. 



4 DISCUSSION AND CONCLUSIONS 

We have found that the SO period of ~ 300 d, which onset in 2005 
was found by R08, has persisted in the X-ray and radio light curves 
of Cyg X-l at least until the switch to the soft state of Cyg X-l in 
2010 June. The statistical significance of this periodicity, at P sup ^ 
310-320 d, corresponds to a chance probability of < 3 x 10~ 5 . 

In addition to the X-ray and radio light curves, we find the 
SO modulation to also be very significantly present in the spec- 
tral hardness, or, equivalently, a local spectral index, in soft X- 
rays, E<10 keV. Also, the depth of the soft X-ray orbital modu- 
lation (caused by absorption in the stellar wind) is found to be a 
strong function of the SO phase, with the orbital modulation being 
much stronger around the minimum flux of the SO cycle than at its 
minimum. Other known cases of correlations between the orbital 
modulation and super orbital phase are those in the low-mass X-ray 
binary 4U 1820-303 I Zdziarski et al. 2007bJ) and in some X-ray 
pulsars (Raichur&Paul 20081) . In the latter case, this phenomenon 
appears also be due to a changing viewing angle of a precessing 
system. 

On the other hand, thanks to the SwiftfBKI data, we find a 
new feature of the SO modulation of Cyg X-l. Namely, the average 
spectral slope at E > 20 keV is completely constant with the SO 
phase. As discussed, e.g., in Z02, the long-term variability of Cyg 
X-l in the hard state is driven mostly by variable irradiation of 
the hot, Compton-scattering plasma, by some soft photons, causing 
changes of the spectral slope and related to changing accretion rate. 
This type of variability is absent when averaged over the SO period. 
This provides strong evidence that the SO modulation is not driven 
by a variable accretion rate. The SO modulation consists of changes 
of the amplitude of the intrinsic spectrum, and of enhanced soft X- 
ray absorption around the minimum flux of the SO cycle. These 
features are fully compatible with the model of a bulge around the 
centre of a precessing disc, proposed by PZI08. 

We note that the precession mechanism should be compat- 
ible with the change of its period by a factor of 2. This ap- 
pears easier to achieve with radiation driven precession than with 
that caused purely by the tid al forces exerted by the donor (e.g., 
IWiiers & Pringlel l999; ODOh lCaproni et alj2006h . OD01 find that 
Cyg X- 1 lies in a region of marginal stability of radiation induced 
precession, which possibly can explain the observed large change 
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of P gup as a manifestation of a chaotic behaviour. Still, if this is the 
case, the rather coherent ~150-d and ~300-d SO modulations over 
time scales of years are difficult to understand. 

An alternative is that the SO modulation is not chaotic but sta- 
ble with the period of ~300 d, but there are changes in the nodal 
precession mode. The observed period of ^ 300 d would then cor- 
respond to the main, retrograde mode with m = bending nodes 
(OD01). The ~150-d observed period would then correspond to the 
prograde m = 1 mode (OD01), in which the actual precession pe- 
riod would be still ~300 d, but the disc bending gives rise to two 
maxima/minima of the observed flux per period. An argument in 
favour of this interpretation is that the precession during the ~150- 
d modulation was found prograde by PZI08. In this model, the ab- 
sorption periodically varying with a half of the precession period 
needs to be due to some medium around the bent disc. 

We have found that the ASM count rates averaged over the 
hard states characterized by P sup ^ 150 d and P<, up = 300 d are 
almost identical. The bolometric fluxes in these two epoch are 
also found to be almost identical (Zdziarski et al ] l201lh . Thus, the 
change of the SO period is not associated with a significant change 
of the average luminosity (most likely tied to the accretion rate), 
which appears to present a difficulty for the model of radiation- 
driven precession. 

We also study the orbital modulation of the X-ray light curves. 
We find that the depth of this modulation in the soft X-rays is signif- 
icantly higher during 2006-2010 than before. This also manifests 
itself in the frequency of the X-ray dips, much higher in this epoch 
than before. This appears to indicate that the stellar wind, for some 
reasons, has become more anisotropic than before. Likely, the bulge 
at the disc centre may have now a higher optical depth than before. 

We find statistically significant orbital modulation of the flux 
(but not of the spectral slope) in the Swift BAT light curve. The 
cause of the modulation appears to be Compton scattering of pho- 
tons away from the line of sight (L06). As the optical depth of the 
wind is higher at the superior conjunction than at the inferior one, 
this scattering causes changes of the observed hard X-ray flux by 
about ±2 per cent. 

Finally, we find no correlation of the X-ray and radio prop- 
erties of the source with the reported detections in high energy 
y-rays. The TeV detection happened during a relatively high X- 
ray flux (see also lMalzac et alj|2008l) . whereas the GeV ones took 
place when both the X-ray and radio fluxes were low. W e note that 
Ferm i LAT detected no GeV emission from Cyg X-l (Hill et al. 
20 id) during the first two fla res reported by AGILE ( Sabati ni et al. 
2010al ; lBulgarelli et alfcOldh . 
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Equation dA3t can be readily solved by Newton's method, in which 
case the derivative is also needed, 



f(g) = a g Rlna + b s lnb. 



(A5) 



The uncertainty of T is given by Ar = AR/\dR/dT\, where AR 
is the uncertainty of R. We find, 



a g (b g - l)lna + ^(l -a g )lnb 

(1 - as) 2 ' 
(In a - In b) In b 

2 In a 



T*2; 

r = 2. 



(A6) 



If R is given by the ratio of photon, rather than energy, fluxes, 
then 2 - T is replaced by 1 - T in equation JAU and g = 1 - T. The 
special cases of T = 2 are now for T = 1 . 



APPENDIX A: CALCULATION OF THE EQUIVALENT 
SPECTRAL INDICES 

Often, hardness of spectra is described by the ratio of count rates or 
fluxes, R, in two channels. However, even if energy or photon fluxes 
are used, this value depends on the width of the channels, and is just 
specific to a given instrument. On the other hand, a hardness ratio 
can be uniquely converted to the equivalent spectral index, T, cor- 
responding to a power law spectrum with the same ratio of fluxes 
as the actual data. Thus, this gives a local spectral slope, which can 
usually be interpreted physically, unlike the hardness ratio, which 
is an abstract and instrument-dependent quantity. Note that this in- 
dex is a unique and monotonous function of the flux hardness ra- 
tio. This index is analogous to colours used in stellar astrophysics, 
which also describe a spectrum in an instrument-independent way. 
This technique was introduced by Z02, who, however, gave no de- 
tails of the calculation. We present them here. 

The condition of the ratio of the energy fluxes in a power law 
spectrum being equal to that in the observed spectrum, R, can be 
written as, 

£f r - £f r 
rf-r p 2-r ' r £ 2; 

£ 2 " L i (Al) 



In E 4 /E 



lnE 2 /Ei ' 



r = 2, 



where the channel boundary energies are (E\, E 2 ), (£3, £4). For 
adjacent channels with increasing energy, £2 = £3 (assumed in 
this work), which, however, is not required in this formalism. We 
define, 

a=§-<l, b=§->l, g = 2-T. (A2) 
In order to find T(R), we need to solve, 

f(g) = b s + Ra g - 1 - R = 0, g * 0; (A3) 

R=~, 8 = 0. (A4) 
In a 
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